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ABSTRACT 


We studied the unique kinematic properties in massive filament G352.63-1.07 at 10?- 
AU spatial scale with the dense molecular tracers observed with the Atacama Large 
Millimeter/submillimeter Array (ALMA). We find the central massive core M1 (12M) 
being separated from the surrounding filament with a velocity difference of v — Usys = 
—2kms'! and a transverse separation within 3 arcsec. Meanwhile, as shown in multiple 
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dense-gas tracers, M1 has a spatial extension closely aligned with the main filament and 
is connected to the filament towards its both ends. M1 thus represents a very beginning 
state for a massive young star-forming core escaping from the parental filament, within 
a time scale of ~ 4000 years. Based on its kinetic energy (3.5 x 1044 erg), the core escape 
is unlikely solely due to the original filament motion or magnetic field, but requires more 
energetic events such as a rapid intense anisotropic collapse. The released energy also 
seems to noticeably increase the environmental turbulence. This may help the filament 


to become stabilized again. 


Keywords: Star formation (1569); Young stellar objects (1834); Dense interstellar 
clouds(371); Interstellar filaments(842); Gravitational collapse(662) 


1. INTRODUCTION 


Star-forming clouds can cause several promi- 
nent dynamical consequences. One is the escap- 
ing and spreading of young stars towards the 
surrounding space, which has caught increas- 
ing attention during the recent years (e.g. To- 
bin et al. 2009; Kraus et al. 2017; Herczeg et al. 
2019; Cantat-Gaudin et al. 2019; Swiggum et al. 
2021). The stars tend to be gradually drift- 
ing over the entire cloud, forming a similar but 
more extended spatial distribution (Grofschedl 
et al. 2018; Jerabkova et al. 2019; Kuhn et al. 
2020; Ward et al. 2020; Gupta & Chen 2022). 
This provides a unique proxy to inspect the 
cloud morphology and mass distribution. Fur- 
thermore, based on the recent Gaia-based mea- 
surement of the trigonometric parallax and stel- 
lar proper motion, one can explore the spatial- 
velocity structure with an unprecedented accu- 
racy and a historical perspective (e.g. Szegedi- 
Elek et al. 2019; Roccatagliata et al. 2020; Galli 
et al. 2020; Krolikowski et al. 2021; Kounkel 
et al. 2022; Tu et al. 2022; Ha et al. 2022; Dhar- 
mawardena et al. 2022; Duan et al. 2023). These 
works have largely improved our understand- 
ing of the dynamical condition and evolutionary 
trend of the star-forming clouds. 

As one critical but incomplete aspect of this 
study, the initial condition of the escaping 
young stars is still not widely inspected. With 
the exception of some extremely high-velocity 


stars (10 to 10? kms !) ejected from n-body 
(n > 3) interactions (Ducourant et al. 2017; 
Fernandes et al. 2019; Bally et al. 2020; Rivera- 
Ortiz et al. 2021), the stellar motions could 
also be affected by their dense-gas structures. 
Young stars could directly inherit the turbulent 
field of the gas structure (Ha et al. 2021; Kro- 
likowski et al. 2021; Quintana & Wright 2022; 
Gupta & Chen 2022), and obtain further accel- 
eration due to the gravitational instability or 
other interactions among the parental filaments 
or clumps (Stutz & Gould 2016; Getman et al. 
2019; Zamora-Avilés et al. 2019; Kim et al. 2019; 
Sharma et al. 2020; Álvarez-Gutiérrez et al. 
2021). Stutz & Gould (2016) proposed one par- 
ticular scenario where the filament in Orion A 
could have an oscillation that drives the young 
stars towards transverse directions, eventually 
forming a broadened stellar distribution in par- 
allel with the filament. But in all these studies, 
the young stars are either totally separated from 
the filaments or have dissipated the surrounding 
gas (e.g. Jerabkova et al. 2019; Gupta & Chen 
2022). So we still need to investigate the initial 
dynamical condition when the YSOs are leaving 
the parental structures. 

G352.63-1.07 is a massive young star form- 
ing region at a precisely measured distance of 
D = 690 pc (Chen et al. 2021, Chen21 here- 
after). It contains several massive cores aligned 
on a dense filament. As shown in Chen21, the 
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filament has actively ongoing dynamical fea- 
tures, including prominent outflows and mass 
transfer flows, which induced bright shock emis- 
sions. However, due to the limited molecular 
tracers and sensitivities, the gas motions among 
the cores and filaments are still undetermined. 
In this work, we explored the gas motions us- 
ing optically thin lines, which demonstrate the 
unique tendency of the escaping massive core 
M1. The observed data is described in Section 
2. The global dense-gas and velocity distribu- 
tion is reported in Section 3. The dynamical 
origin and the time-energy properties of the gas 
motion are discussed in Section 4. A summary 
is given in Section 5. 


2. OBSERVATION 


G352.63-1.07 (G352 hereafter) is ob- 
served with the Atacama Large Millime- 
ter/submillimeter Array (ALMA) in the 
ATOMS survey (Liu et al. 2020, 2022). The 
observation was performed with both Compact 
7-m Array (ACA) and the 12-m array (C43- 
2 or C43-3 configurations) in Band 3. The 
restored data cube has an average beam size 
of Omaj X Omin = 2”.5 x 21 (PA = er 
The antenna baselines can cover the extended 
structures with a scale up to ~ 100”. The 
channel width and noise level varies with dif- 
ferent spectral windows. The H!?CO*, CCH, 
and HPCN (1-0) lines have Ava, = 0.2 km s7}, 
Orms = 5 mJy beam! (0.15 K). The HCO* (1-0) 
has Augg = 0.1 kms !, o4, = 9 mJy beam! 
(0.25 K). The CHOH (211 — lio) lines have 
Ucan = l.4kms !, Orms = 2mJy beam"! (0.08 
K). The more detailed observing conditions are 
presented in Liu et al. (2020). 


3. RESULTS 


3.1. Filament and cores 


Figure 1a shows the 3 mm continuum emis- 
sion and the integrated HCO* (1-0) emission 
(moment 0). Compared to the chemically fresh 
molecules CCH and HC3N (Chen21, Figure 3 


therein), the H^CO* emission is more concen- 
trated towards the inner region of the filament 
structures. It shows an S-shaped compact fila- 
ment. The six major cores (Chen 21) are clearly 
resolved therein. The line intensity decreases to- 
wards M4 and M5, which could be due to the gas 
dissipation caused by their outflows (Appendix 
B). 

Figure 1b shows the H^CO* emission over- 
laid on the Spitzer/IRAC three-band color im- 
age (GLIMPSE survey, Benjamin et al. 2003; 
Churchwell et al. 2009). Among the dense cores, 
only M1 has a compact IR source. M2 to M5 
are all absent of IR point sources. The other 
two IR. sources (M3b and M6) are located on 
the western side of the filament. The IRAC 
color magnitudes for M1 are [3.6]-[4.5]=2.1 and 
[5.8]-[8.0]=1.0, which is similar to the color of 
Class-0 YSOs (Megeath et al. 2012). M1 is also 
detected in 6.7 GHz methanol maser but un- 
detected in 6-8 GHz radio continuum emission 
above the noise level of 3-8 mJy beam! (Walsh 
et al. 1998). This indicates M1 to have a deeply 
embedded young massive star. It has not yet 
caused significant ionization to the surrounding 
medium. Figure 1b also shows the CH3OH (2-1) 
emission. It has a compact morphology concen- 
trated at M1. 

The molecular spectra at the three inner cores 
are shown in Figure 1c. The H'8CO* (1- 
0) exhibits a noticeable double-peak profile at 
M1. The two components have a separation 
of Ubiue — Ura = —1.8kms ^. M2 and M3 
both have a single-velocity component within 
Ug; = £0.3kms~', which are close to the red- 
shift component at M1. The HCO* and CH3OH 
line profiles are less resolved owing to high opti- 
cal depth and low spectral resolution, respec- 
tively. But we can still see these two lines 
inclined to the blueshift side. We also exam- 
ined the other two dense-gas tracers CCH and 
HP?CN (1-0) (Figure 3). They have similar ve- 
locity components to the H?CO* (1-0). The 
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molecular lines all indicates prominent blueshift 
component at M1. 


3.2. Physical parameters and gravitational 

instability 

The physical parameters of the cores are pre- 
sented in Table 1. The details of the calcula- 
tion are described in Appendix A. All the cores 
have comparable spatial sizes are nearly uni- 
formly distributed along the filament, with an 
average spatial interval of Asgores = (10 + 2)", 
or (7+1.5) x 10? AU. 

From the observed parameters, we estimated 
the stability of the cores from the critical mass 
(Bertoldi & McKee 1992; McKee & Ostriker 
2007; Li et al. 2013): 


Merit = Mge + Mo 
= z (s. $ 4) "EL 


wherein Grot is the effective velocity disper- 
sion (see Appendix A), Mpg and Ma represent 
the upper-limit masses to be supported by the 
turbulence and magnetic pressures, respectively. 
Ma depends on magnetic field strength B and 
mass density po = mino, which are included 
in Alfvén velocity of v4 = B/\/4mpo. For the B 
value, we referred to the B-field measurement 
in other similar cold dense filaments (e.g. Liu 
et al. 2018; Ching et al. 2018). They are mea- 
sured to have a bulk distribution of B — 0.5 
to 1.0 mG. Adopting this range, we can derive 
Ma = 0.6 to 1.2M5, which only has a minor 
contribution to Merit- As shown in Table 1, all 
the cores have Merit > Meore, which indicates a 
subcritical state if the turbulence can support 
them against the self-gravity. 

The filament stability can be estimated from 
the critical line-mass density (Ostriker 1964; Ar- 
zoumanian et al. 2013) 


M 207 4 
— E a. 2 
( l ).. G l ) 


From the average line width of 
Avgia(H®CO*) = 2.0kms !, we can derive 
Otot = 0.9 kms! and (M/l) ait ~ 390 Mo pc. In 
comparison, the observed total mass (48 Mọ) 
and length (45") yield (M/l)as, ~ 300 Mo pc !. 
Considering the projection effect, the actual 
value could be even smaller. Therefore the 
property of (M/l)obs < (M/l)a& suggests a 
subcritical state also for the entire filament. 
The turbulent condition of the filament and 
cores should have a close interplay with other 
dynamical features including the core collapse 
and escaping motion, as discussed bellow. 


4. DISCUSSION 
4.1. Resolving the velocity components 


As shown above, the filament has an overall 
smooth and compact morphology. The transfer 
flows (Chen21) and outflows seem to only have 
a limited influence to the main filament. As 
seen in optically thin lines, the blueshift motion 
of M1 is still the most conspicuous feature over 
the entire structure. Its dynamical origin should 
be further inspected. 

Figure 2a shows the H'?CO* (1-0) emission 
in four different velocity channels. The ALMA 
3 mm and SMA 1 mm (Chen21) dust contin- 
uum emissions are also plotted in sub-panels for 
comparison. Figure 2b presents the position- 
velocity (PV) plot of the H?CO* along the fil- 
ament. It shows that the blue component is 
mainly in the velocity range from -4 to -1 km s7! 
and has a spatial extension from offset— — 106" to 
+6”. In order inspect the gas motion, we plot 
the emission region in two velocity intervals, de- 
noted as low- and high-velocity components, or 
LVC and HVC, respectively. we plot their emis- 
sion regions in Figure 2a. One can see that the 
LVC is peaked at M1 and has a weak elonga- 
tion towards M2, while the HVC is more con- 
fined around M1. Their morphologies are both 
largely different from the outflow lobes in Fig- 
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ure 6a. They should both trace the dense-core 
motion instead of outflow. 

The filament component is mainly in the ve- 
locity range of (-1.5,+1.5) kms~!. Its emission 
region (false-color image in Figure 2a) almost 
traces the entire gas structure, and shows a no- 
ticeable gap at M1 center. The redshift wing 
component (1.5 to 2.2 kms !, yellow contours) 
includes two separated patches, with one ex- 
tending from M4 and M5 and another around 
M3. They could trace the denser gas affected by 
the outflow or transfer flow. Around M1, there 
are no evident redshift feature at v > 1 kms. 
The dust continuum emissions (Figure 2a, two 
right panels) show an overall similar spatial ex- 
tent with the H'3CO?* filament, but have a com- 
pact intensity peak rightly at M1, which is ~ 3 
times more intense than the remaining filament. 

The spatial correlation of the velocity compo- 
nents can also be inspected from their intensity 
profiles as shown in Figure 2b (lower panel). 
The dust continuum profiles (dashed and dot- 
ted lines) rightly follow the blueshift compo- 
nent (blue solid line). They all have a shoulder- 
like decreasing trend from M1 to M2. This co- 
herency suggests that the dense core M1 should 
be mainly associated with the blueshift gas. 

The intensity gap on the filament is more 
clearly seen in the intensity profile (red line in 
Figure 2b, lower panel). It rightly coincides 
with the M1 center of both the dust contin- 
uum and the blueshift H!?CO* profiles. For 
the HP?CO* emission, if considering the total 
intensity profile (black solid line), the blueshift 
component would nicely fill the gap, making the 
entire profile much more flattened. 

The CCH emission regions of the two com- 
ponents are shown in Figure 3a. The PV dia- 
grams of CCH and HP?CN are shown in Figure 
3b. The lower intensity around M1 along the 
filament can also be seen in CCH (Figure 3a). 
Although its blue component is more extended 
than that of H?CO*, it is still mainly concen- 


trated at the central cavity of the filament. The 
H!?CN emission is more concentrated at M1, so 
the central intensity decline is not evident. But 
from its PV diagram, one can still see a compa- 
rable velocity separation of v — vs, = —2 km s7! 
between the blueshift component and the main 
filament. And it is also closely overlapped with 
the H?CO* and CCH spectra for the blueshift 
peak (Figure 6c). The three molecular lines thus 
consistently suggest the blueshift component to 
be the dominant one at M1. In other words, the 
dense core is having a bulk motion relative to 
the rest part of the filament. 

Based on their spatial and velocity features, 
M1 and the main filament could have two possi- 
ble configurations. M1 can be either leaving the 
filament on the front side or moving towards it 
from behind. From the PV plot (Figure 2b), one 
can see the blue component connected to the 
main filament both towards north and south, 
around offset —4-5" and —15", respectively. The 
PV diagrams of the CCH and HPCN (1-0) lines 
(Figure 3, right column) show the similar con- 
nection feature between M1 and the main fil- 
ament in the two directions. M1 should thus 
be originally formed in the filament, and have 
obtained the blueshift motion only during the 
recent time. 


4.2. Driving force of the blue component 


The multiple velocity components in G352 
are comparable to other filaments with promi- 
nent kinematical features. As shown in previ- 
ous studies, undisturbed linear filaments tend 
to have moderate fluctuation of |v — vy,| < 
l.O0kms ! around their internal cores (e.g. 
Punanova et al. 2018; Bhadari et al. 2020). 
More complicated fiber-composed filaments also 
have similar velocity variation within lkms'! 
(Hacar et al. 2017; Clarke et al. 2018). Stronger 
kinematical features of several kms~! are seen 
in filaments with collapse (Henshaw et al. 2014; 
Chen et al. 2019; Ren et al. 2021; Li et al. 2022; 
Cao et al. 2022) or interactions (Shimajiri et al. 
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2019; Anathpindika & Francesco 2021). The ve- 
locity fluctuations in G352 resembles those most 
active ones. The blueshift motion (—2 kms~') 
and its spatial scale (5" or 3500 AU) lead to a 
gradient of ~ 250kms~'pe~!, which resembles 
those most intensely collapsing filaments (e.g. 
Peretto et al. 2013; Montillaud et al. 2019; Hu 
et al. 2021; Beuther et al. 2021; Cao et al. 2022). 

In collapsing or fiber-composed filaments, the 
global velocity field is often entangled with the 
individual core motions. It is therefore uncer- 
tain if the cores are co-moving with the fil- 
ament or already separated. In comparison, 
G352 presents an example of clearly separated 
velocity components. In fact, the filament is 
also inclined to the blueshift side around M1 
towards -0.5 kms ! (Figure 2b). The filament 
could thus have a co-moving tendency with M1. 
In particular, as shown in Figure 2a, the HVC 
(v = (—4, —3) kms!) is narrowly confined be- 
tween the northern and southern segments of 
the entire filament. This increases the evidence 
that the core and filament motions are closely 
related. 

The mass transfer flows (Chen21) provide a 
viable mechanism to initiate the filament col- 
lapse. Since the flows are observed from one- 
sided molecular line wings instead of strong in- 
fall signatures, they would provide a moderate 
mass accumulation to the inner filament (be- 
tween M2 and M3). If the mass assembly once 
exceed the threshold of (M/l)cit or Merit, it 
would possibly induce a major collapse. One 
can see that the transfer flows exhibit no evident 
mass assembly at their arrival points on the in- 
ner filament (Figure 2). This also indicates that 
the flows should have a further propagation to- 
wards M1. 


4.3. Energy scales of the gas components 


The filament collapse and core escape can also 
be examined from their energy scales. From 
the core mass and its escaping velocity, the ki- 
netic energy of M1 is estimated to be Ek = 


Meore(Ublue — Vrea)” /2 œ 3.5 x 10“ erg. This 
would represent a lower limit due to the projec- 
tion effect. In comparison, the turbulent energy 
of the entire filament is Ei = mQc2,/2 = 
(2.5 3- 0.5) x 10“ erg. As for the magnetic en- 
ergy, from the speculated B-range (Section 3.2), 
we can estimate Eg = m,4,U2/2 = 2 to 8 x 10% 
erg. Ej, and Eg values characterize the en- 
tire filament. The fraction to reach M1 could 
be even lower thus cannot be solely responsible 
for the core motion. 

The collapsing energy can be estimated as 
AEp co = Gm? [(1/r1) — (1/ro)], which involves 
the initial and final radii (ro and rı) of the col- 
lapsing mass "my. Assuming that the major 
collapse is taking place between M2 and M3, 
we can adopt ro = AsScores. Together with 
the current radius and mass of M1, we derived 
AE, cou c 6 x 10“ erg, which could be sufficient 
to accelerate M1. 

As seen in PV diagram (Figure 2b), the in- 
ner region has a typical line width of Av > 
l.8kms^! for both the blueshift component 
and the main filament. It rapidly declines to- 
wards the outer part, reaching a trans-sonic 
level of Av ~ 0.7 kms ! (o, = 0.4kms 1) 
towards offset=+20”. The low-Av areas over 
the filament could represent the fraction less af- 
fected by the collapse. Adopting the average 
value (0.9 kms!) in Equation 3, we can derive 
(M/l) ait & 150 Mo pe, which is much smaller 
than (M/l)obs and should reflect the supercriti- 
cal condition before the major collapse. 

Figure 4 shows the velocity distributions over 
the core and filament in all three lines. The low 
turbulence areas with Av down to 0.6 kms"! 
can be seen in all three species. The CCH and 
H'3CO* emission regions both exhibit the low 
values down to Av ~ 0.7 kms^!, while the 
higher values mainly appear between M1 and 
M3. The H!3CN line is more concentrated on 
dense cores, which could be due to its chemi- 
cal bias towards the protostellar stage. Along 
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the filament, the H'°CN line width variation 
is mainly between 1.2 and 1.4 kms !, which 
is comparable to average values of CCH and 
H?CO* lines. 


4.4. Time scales of the gas motion 


Figure 5a shows the ridge lines of the main 
filament and the blueshift component (M1 and 
southern extension). The two components have 
a transverse offset of As... = 0 to 2.3 arc- 
sec along the filament. We plot the As, and 
lUbiue — Urea| profiles along the filament in Figure 
5b. One can see that the two quantities have a 
coherent variation trend. At M1 and the sec- 
ond velocity peak near M2 (offset——10"), the 
velocity and transverse offsets both reach a lo- 
cal maximum. 

The transverse separation as a function of 
time can be estimated as Asec = tesc(Ublue — 
Urea) tan 0, wherein 0 is the inclination angle be- 
tween the gas motion and sightline. Adopt- 
ing 0 = (45 + 10)°, we can derive a timescale 
of tec c (4+ 2) x 10? years for the blueshift 
motion. It is notably close to the outflow age 
of tout = 3.6 x 10° years (Appendix A). The 
two values should together characterize the star- 
forming age in G352. 

If the core escape started during the filament 
collapse, one would expect an overall collapsing 
time sufficiently longer than tese. From the re- 
cent semi-analytical work (Clarke & Whitworth 
2015), the filament could have a collapsing time 
of 

teoll = (0.49 + 0.26A9)(Gpo) ^, (3) 


wherein A is the aspect ratio of the filament. 
Assuming that the major collapse took place 
in the inner region between M2 and M3, from 
the filament length (| ~ 22") and width (w ~ 
5"), we can derive teon c 6 x 10* years. It 
could roughly represent the dense-core forma- 
tion time. For M1, it then implies a mass 
accretion rate of Mace = Mee) toon ~ 1.8 x 
1074 Mo yr-* ^, which is comparable to the rate 


of the transfer flows (chen21). The collapse to- 
wards center could be dominated by the two col- 
liding flows along the filament. The flows would 
rapidly increase the central density and develop 
an anisotropic instability. The highly com- 
pressed gas would then move towards the trans- 
verse directions, initiating the observed motion. 
A schematic view of the collapse and core accel- 
eration is shown in Figure 5c. 

The collapsing time could be compared with 
the turbulence-increasing time of the filament. 
From the turbulence level of the inner filament 
(on ~ 0.8kms7!), we can derive a dynamical 
timescale of tayn = AScores/Ont ~ 8 x 104 years, 
which is also similar to teo. As shown in Fig- 
ure 4, the line-width distributions of H^CO* 
and CCH are both steeply increased between 
M1 and MB with a scale of AV > 1.0kms-!. 
This also provides evidence that the turbulence 
was increased by the recent dynamical evolution 
thus have not yet reached a uniform state over 
the filament. 

Although the filament collapse could have suf- 
ficient energy and time to initiate the core mo- 
tion, other dynamical origins for the blue com- 
ponent are still not fully excluded. In partic- 
ular, it could also be a fiber-like component. 
The fibers in massive filaments are observed on 
larger spatial scales (> 0.2 pc) in a few re- 
cent studies (Shimajiri et al. 2019; Cao et al. 
2022). Those fibers are only loosely aligned or 
intertwined along the filament axis. They have 
transverse separations up to several 0.1 pc, and 
do not have two-ends connections like in the 
case of M1. The fibers could be increasingly 
separated with the filament age. For G352, the 
separation would be Ass > (Uplue — Ured tco & 
2.5 x 104 AU (36”), which is unreasonably large 
compared to the observed spatial scales. To 
maintain the observed Asesc, the fibers should 
have originally much smaller velocity difference 
(« 2kms 7), and still need a recent collapse to 
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provide the kinetic energy and reach the escap- 
ing velocity. 


5. SUMMARY 


Based on the velocity distribution of H'?^CO* 
and other dense-gas tracers, we found that the 
most massive core M1 (12 Mọ) is connected 
to the main filament, but has a blueshift of 
Vcore — Usys X —2 km s7}, The core motion has 
a kinetic energy of 3.5 x 104 erg and timescale 
of only ~ 4000 years. As a unique example, M1 
reveals the very beginning of a massive YSO es- 
caping from its birth site. It confirms the pos- 
sibility that YSO can leave the filament along 
the transverse direction. 

Among the available mechanisms, only the fil- 
ament collapse could provide enough energy to 
initiate the core escape. A recent collapse is also 
needed to generate the massive star therein. It 
can be also responsible for the drastic velocity 
and turbulence variations over the filament. As 
the core motion is taking place, a central gap is 
formed rightly at its location on the remaining 
part of the filament. 

According to the energy and time compar- 
isons, the process of core escape could have 
three major steps: (i) the filament had a rel- 
atively low turbulence, with moderate contrac- 
tion and fragmentation into the dense cores; (ii) 
because of the transfer flows, the inner region 
became more massive and unstable, and even- 
tually initiated a major collapse; (iii) M1 was 
strongly compressed by the collapsing gas to 
start the escaping motion, while the rest energy 
could be transferred to the filament to increase 
the turbulence. 

Based on the result in G352, one can exam- 
ine other filaments with steep velocity gradients 
to see if similar core motions are taking place. 
A larger sample will help evaluate the univer- 
sality of such collapse-induced YSO escape, and 
estimate its contribution to the YSO dispersion 
over molecular clouds. 
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Table 1. The physical properties of the cores. 


Parameters M1 M2 M2b M3 M4 M5 
(Observed) 

Ug. (HCO*) (kms!) -2.5 -0.8 -0.6 +0.1 +1.0 +1.3 
Tex (HOO*) (K) 37 22 22 22 22 20 
T, (H'3CO*) (K) 4.0 3.0 5.0 2.1 3.7 2.6 
Av (kms™!}) 2.5 2.3 1.8 2.2 1.8 17 
Radius (arcsec) 5 4 6 5 4 4 
(Derived) 

Niot (1073 cm7?) * 26403 22402 2340.3 24402 22402 22402 
Otot (kms 1) 1.1 0.9 0.8 1.0 0.7 0.7 
Mass (Mo) 12 +2 642 11 E2 8+2 T2 T2 
Merit (Mo) 15 13 12 12 11 11 


a. Average radius deconvolved with the beam size. 

b. Nis. is derived from H!?^CO* intensity using Equation A4. For M4 and M5, the 
H!?CO* emissions are noticeably dissipated by the outflows. We assumed them to 
have same Ntot with M3 based on the fact that these cores also have comparable 3 


mm continuum emissions. 
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Figure 1. (a) The velocity-integrated intensity of the H'’CO* (1-0) line (false-color) and 3 mm continuum. 
The contour levels are 40;ms (1.6 mJy beam~') to 840;ms (peak) in step of 160;ms. (b) The H'?CO* and 
CH3OH emission regions around the main filament, overlaid on the IRAC-RGB image (3.6, 4.5, and 8.0 
um bands). The H1?CO* contours are 15% to 9096 in 15%-step of the peak intensity (8.5 K kms~‘). The 
CH3OH contours are 1096 to 9096 in 20%-step of the peak intensity (18.7 K kms"). The dashed circles 
labels the area of each core. (c) The spectra at the selected core centers. The blue and red-shaded areas 
indicate the velocity ranges of the two velocity components, respectively. The vertical dashed line denotes 
the division between the blue and main-filament components. 
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Figure 2. (a) Contours: emission region of the H!?CO* (1-0) line in three velocity intervals. The background 
image is the HI^CO* emission in (-1.5,+1.5) km s^! (red-center component). For the blue-wing component, 
the contour levels are 4, 6, and 8 times of the rms level (0.3 K kms^!). For the other components, the 
contour levels are 10% to 90% in 20%-step of the peak intensity, which is 15, 11, and 3 K kms^! for the 
blue-center, red-center, and red-wing, respectively. The green arrows label the mass transfer flow directions 
onto the main filament (Chen21). The green circles labels the possible arrival points of the transfer flows 
onto the filament. (b) PV plot and intensity profile along the major axis of the main filament. The sampling 
direction is labelled in dashed line in panel (a). The vertical dashed lines denote the projected offset of 
the three dense cores on the sampling direction. The horizontal dotted line represents the average systemic 
velocity of the filament. 
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Figure 3. Left column: The emission regions of the blue- (contours) and redshifted (false-color) components 
in CCS and H!8CN (1-0) lines. The contour levels are 20% to 90% of the peak intensity. Right column: The 
PV diagrams of the two molecular lines. the H?^CN emission at v, = 2 to 4 kms~! is from anther HFC of 
F —1-—1. The HFCs of H?CN are separated for 6-7 kms~! and would not blend with each other. The 
CCH (1-0) emission shows an additional small blueshift wings around M3, which should correspond to the 
transfer flow onto the filament (Chen21). This feature is not seen in other lines probably because of their 
lower optical depths. 
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Figure 4. The peak velocity and line width distributions of the two components as observed in three 
molecular lines. The velocity ranges to measure the blueshift and main-filament components are (—4, —1.8) 
and (—1.5,2) kms~?, respectively. left two columns: The radial velocity distributions of the two components. 


right two columns: The line-width distributions of the two components. 
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Figure 5. (a) The ridge lines of the two velocity components, blue component in blue line and dots, red- 
center component (main filament) in yellow line and dots. The dot size is proportional to the integrated 
intensity. The arrows denote the projected moving direction of each component expected on the sky-plane. 
(b) The transverse offset and velocity difference between the ridge lines of the two components along the 
filament axis. In the lower panel, the blue and red dots represent the peak radial velocity profile of the two 
components, respectively. The gray dots represent |vpiue — vrea]. (c) A schematic view of the gas motion on 
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APPENDIX 


A. DERIVING THE DENSE CORE PROPERTIES 


The total column density of a molecular species is estimated from the optically thin lines (Caselli 
et al. 2002; Henshaw et al. 2014) as 


N u ST lix gu 1 
LEE E gTa a] 
" 1 
1 — exp (-hv/kpTa) 
M Du 
gi exp (—El/kpTex)’ 


wherein A, = f Tydv is the total intensity of the line profile, A is the Einstein coefficient for 
spontaneous decay, gu and gı are the statistical weights for the upper and lower states, respectively. 
J(T) = To/|exp(To/T) — 1] is the Planck-corrected brightness temperature with the reference value 
of To = hv/k. Qrot(Tex) is the partition function, A and v are the wavelength and frequency of 
the line transition. Ty, = 2.73 K is the cosmic background temperature. kg and h are the Boltz- 
mann and Planck constants, respectively. The total gas column density is N(H2) = Mtot.moi/Xmol; 
wherein Xma is the molecular abundance. For the H'?CO* abundance, Gerner et al. (2014) measured 
Xga(HPCO*) = 0.9 to 1.5 x 107? from a large sample of dense cores from IRDC to hot-core stage. 
We adopt X,,4(H?CO*) = (1.2 + 0.3) x 107° for the calculation of G352. 

We also tried to use CCH and H!3CN (1-0) lines to calculate Miot- They provide comparable values 
on 10% cm"? scale. The optical depth of the molecular lines are examined using (Bell et al. 2014): 


e Ai GuN tot Eup hvij 
dm Srv? AvQrot PSP (- a) E (=) i T ` (A2) 
In calculation, we also adopted the average abundances in the young high-mass protostellar ob- 
jects (HMPOs) in Gerner et al. (2014), which are X4,4(CCH) = 5.0 x 1078 and X44(HI35CN) = 
Xma(HCN)/80 = 4.7 x 107". From the physical condition at M1 of No, = 2.6 x 10? cm? and 
Tu = 37 K, we estimated 7(CCH) = 1.1, zr(H900*) = 0.12, and 7(H CN) = 0.06. 

From N (H5) distribution, we can estimate the core mass to be Meore = umy f N(H2)dA, wherein 
the integration is made over the core area. The cores are irregular and not fully separated from 
the filament. To prevent complexity, we assumed the cores to have spherical or elliptical shapes. 
M3 and M5 indeed have noticeable elongation thus are assumed to be elliptical, while other cores 
are considered to be spherical. The core region is adjusted to cover the emission region above the 
extended emission over the filament (3 K kms~!). The core areas are plotted in Figure 1a. 

In deriving Noi, we adopt the excitation temperature of HCO* (1-0) lines (Figure 1). Assuming 
optically thick, it can be estimated using 


(A1) 


Testo qu. (A3) 


The calculation leads to Tex ~ 37 K for M1 and 22 + 3 K for other cores. The second value is 
comparable to the dust temperature derived from SED fitting (Chen21), whereas the higher value 
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at M1 should reflect the stellar heating. From CH3CN population diagram, Chen21 derived a much 
higher temperature of Tiot = 165 K at M1. This value should incline to the small amount of shocked 
gas, whereas the major fraction of the filament should still be much cooler due to the absence of 
ionized gas and strong IR emission. Their temperature range is comparable with the values of pre- 
stellar cold dense cores (e.g. Xie et al. 2021), suggesting an early evolutionary stage on average for 
the filament. Moreover, if assuming T4,(H'^CO"*) = 165 K, one can derive N(H2) = 7 x 10? cm? 
and Meore = 30 Mo for M1. This value is unreasonably high as it takes up more than 1/2 of the total 
filament mass. 
The total velocity dispersion of the core is estimated from the observed line width Av as 


"m kpTkin Av? ket kin (A4) 
Vo uma 812 Mma’ 
wherein Mo) is the molecular mass. The thermal part is 
k Ti in 
on — 4/—— (A5) 
HMH 
and the non-thermal part is 
Av? kpTin 
n = m P A 
oe "Ts Mmol (A6) 


In calculation we adopted a kinetic temperature of Tij, = 37 K for M1 and Tkin = 22 K for other 
cores. They lead to cq, = 0.29 to 0.36 kms !, respectively. In comparison, the average line width of 
Av = 1.5km s™! leads to on, = 0.64 kms" !, which is not largely affected by Tin and o. 


B. OUTFLOW PROPERTIES 


The outflow emission is detected from high-velocity line wings in the HCO* (1-0) lines. Its emission 
region and spectra are shown in Figure 6a and 6b, respectively. The blue- and redshift line wings are 
both strongly detected at M4, which shows the velocity ranges of (—20, —6) kms~! for the blue lobe 
and (5,20) kms! for the red one. As seen in Figure 6a, compared with the CO outflow (Chen21), 
the HCO* is concentrated around M4 and M5, suggesting that the outflow should be launched from 
this area. From the spatial extensions of the outflow lobes, we can see at least two bipolar outflows, as 
indicated by the dashed lines. The outflow extensions are not fully overlapped with M4 and M5. It is 
uncertain if this is due to additional driving sources or the time-variation of the outflow morphology. 
One thing for sure is that the outflow intensity is much weaker at M1, and the red lobe is almost 
undetected. If M1 also has a contribution to the blueshift gas, it should be more likely ejected in the 
process of core collapse and escape than due to the protostellar outflow. 

From the HCO* line-wing intensities and outflow emission areas, they are found to have masses 
Mblue = 0.20 Mo and Mreq = 0.24 Mo. And the average radius is Rou = (843)” = (5.542) x 10° AU. 
Assuming an average velocity of Vou = 8kms7!, the outflow age is derived to be tout = Rout/Vout = 
(3.5 +1) x 10° years. 
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Figure 6. (a) The HCO* outflow lobes overlaid on the HCO* excitation temperature map. The line intensity 
is converted to Tex using Equation A3. The contour levels are 15% to 90% of the maximum intensity, which 
is 22.6 and 22.9 K kms^! for the blue and red lobes, respectively. (b)(c) The molecular spectra at selected 
cores. For each spectrum, the horizontal line represents the zero-level base line. The vertical dashed line 
denotes the division between the blue and main-filament components. 
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